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We study the Regge-Wheeler and Zerilli equations (RWE and ZE) at the 'algebraically special 
frequency' Q., where these equations admit an exact solution (elaborated here), generating the SUSY 
relationship between them. The physical significance of the SUSY generator and of the solutions at 
Q. in general is elucidated as follows. The RWE has no (quasinormal or total-transmission) modes 
at all; however, is nonetheless 'special' in that (a) for the outgoing wave into the horizon one 
has a 'miraculous' cancellation of a divergence expected due to the exponential potential tail, and 
(b) the branch-cut discontinuity aX lu = Q, vanishes in the outgoing wave to infinity. Moreover, 
(a) and (b) are related. For the ZE, its only mode is the — inverse — SUSY generator, which is at 
the same time a quasinormal mode and a total-transmission mode propagating to infinity. The 
subtlety of these findings (of general relevance for future study of the equations on or near the 
negative imaginary ti;-axis) may help explain why the situation has sometimes been controversial. 
For finite black-hole rotation, the algebraically special modes are shown to be totally transmitting, 
and the implied singular nature of the Schwarzschild limit is clarified. The analysis draws on a 



recent detailed investigation of SUSY in open systems |math-ph/990903C] 
PACS numbers: 04.30.-w, 04.20. Jb, 04.70.Bw, 11.30. Pb 



I. INTRODUCTION 

The evolution of small perturbations of the Schwarz- 
schild metric of a spherically symmetric black hole, such 
as generated by infalling matter |^ or in the aftermath of 
a stellar collapse Q, is well known Q to be described in 
each angular- momentum sector > 2 by a Klein-Gordon 
wave equation 



[dl+oj"^ -V[x)](t){x,u) =0 



(1.1) 



Here, </) is a scalar field (a combination of the metric- 
function perturbations), while the 'tortoise coordinate' 
is 



X — r -\- ln(r — 1) 



(1.2) 



where r is the circumferential radius (so that 1 < r < oo 
maps to — oo < x < oo), and c — G — 2M — 1 {M is the 
black-hole mass). However, the axial sector is described 
by the Regge-Wheeler equation (RWE), which has the 
potential 



V{x) = (r- 1) 



2(n+l)r- 3 



Clearly, one is now interested in the modes of these 
equations, i.e., eigenfunctions under either outgoing- 
{(f>{x^±oo,uj) ~ e^'^l'^l; however, see below for precise 
definitions) or incoming- wave ((/)(x,w) ~ e-*'^l'^l) bound- 
ary conditions Foremost among these are the quasi- 
normal modes (QNMs, with outgoing boundary condi- 
tions), which constitute the excitation spectrum of the 
black hole. Besides being a matter of principle in decid- 
ing the stability of the Schwarzschild solution, the issue 
takes on exciting significance with the prospect of their 
experimental detection in the coming decade using in- 
struments such as LIGO and VIRGO |||. QNMs also 
play a role in fully nonlinear numerical simulations [Q, 
where they are observed to dominate the radiated signal 
at intermediate times after a violent event such as a stel- 
lar collapse or a black-hole collision (cf. the beginning of 
this section), so that numerical and perturbative studies 
are complementary. Besides QNMs, we will also consider 
total-transmission modes, incoming from the left (right) 
but outgoing to the other side, and denoted as TTMlS 
(TTMrs). The fourth type of mode, a normal mode 
(NM) or bound state with (j){x,Lo) ~ e"!'^'^', is ruled out 
since both V and V are purely repulsive. 



(1.3) 



while in the polar sector, the Zerilli equation (ZE) has 
8n^(n-|-l)r^ + \2n'^r'^ + 18nr + 9 



V{x) = (r- 1)- 



r4(2nr + 3)2 



(1.4) 



where we defined n — i(i?— 1)(^+2). 



It has been noticed ||] that the RWE and ZE have the 
same QNM spectrum, which a priori is by no means obvi- 
ous given the very different forms for the potentials. This 
can be understood because solutions of the two equa- 
tions are related by 'intertwining' ^ or supersymmetry 
(SUSY) ||ll as = [dx + W{x)](j), with the superpotential 



1 



where 



N 



2 r2(2nr + 3) ' 
4n(n+l) „/£ + 2 



(1.5) 
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Of particular interest are now the SUSY generator 

ei(x-) =exp|- / dyW{y) 



(1.6) 



(1.7) 



(cf. Section [IB for the notation) and its counterpart 



= 1 11 1, which generates the inverse transform. By 



the general theory of SUSY, these are solutions of the 
RWE and the ZE respectively, at the same eigenvalue, 
which in this case is yielded by (1.1) as — fl^ for the 
'algebraically special frequency' Q = —iN/2. 

Since the early study of 'algebraically special' black- 
hole perturbations, these have been associated with 
TTMs: Ref. Q, deahng with the Schwarzschild case, 
states (below its last, unnumbered display) that 'The ra- 
diation field (incoming) . . . vanishes.' The analysis of ||l^] 
was generalized to Kerr holes in |l^ , still on the level of 
general relativity, and subsequently Ref. gave the cor- 
responding exact solutions of the RWE, ZE, and Teukol- 
sky equations. Given the fact that |Q consistently calls 
these solutions TTMs, it is remarkable that their solv- 
ing the appropriate wave equations is checked explicitly, 
but not their obeying the corresponding boundary con- 
ditions. 

The situation was confounded when Ref. p5| ] reported 
a QNM il' of the RWE numerically very close to, but 
not at, Q. For £ = 3, however, f2' is given with a nonzero 
real part, while the QNM — fi'* which would then be 
dictated by symmetry is not found; this seems to imply 
that the last 3-4 digits (the order of magnitude of the 
reported Re il') given by are all insignificant, so that 
no conclusions can be drawn about possible coincidence 
of fl and ri' (if the latter does exist at all (l6j). In con- 
trast, Ref. reports TTMrs of the RWE very close 
to for 2 < £ < Q, and suggests that these correspond 
to the special mode; again, if all the given decimals are 
significant, strictly speaking this suggestion is not war- 
ranted. Finally, Ref. |lj] finds that for £ — 2 and for each 
'magnetic quantmn number' —2 < m < 2 the Kerr hole 
has a QNM (the ninth) which tends to = —Ai in the 
Schwarzschild limit a ^ (a is the black-hole rotation). 
However, Ref. [|l8| cautions that these QNMs should dis- 
appear for a = 0, since they cannot coexist at 51 with the 
TTM furnished by the special mode. For each m, this dis- 
appearance is thought to occur by cancellation of the two 
modes (mirror images only for to = 0) on each side of the 
negative imaginary axis (NIA) in the w-plane. Certainly 
this is conceivable, since an analogous cancellation hap- 
pens for the Poschl-Teller potential V{x) = Vcosh^^(.T) 



(sometimes used as a model for gravitational potentials) 

(e-g- 



for ^i-V = ^, §,§,... (e.g., mC 

While the mode situation at fl for the RWE alone thus 
is unclear already, SUSY leaves its relationship to the ZE 
nontrivial at uj"^ = [|l^, spanning the kernel of the 
SUSY operator dx + W. It should be noted that the exist- 
ing literature rarely emphasizes this potential difference 
in spectrum between the two equations: for instance, the 



captions to Figs. 1 and 2 and Table 1 of |15| merely say 
'Schwarzschild QNMs,' although the accompanying text 
makes it clear that in fact the RWE has been studied. 
In view of this state of affairs, it seems imperative to 
start from first principles, using comparatively rigorous 
methods, and study each of the six questions: at f2, do 
the RWE and the ZE have a QNM and/or TTMr and/or 
TTMl? If this investigation can clarify the pitfalls await- 
ing a numerical analysis, it will be of value beyond this 
immediate scope. 

We will settle the issue as follows: at il, the RWE has 
no modes at all, while the ZE has one mode which 
is simultaneously a QNM and a TTMi^. To this end, in 
Section || we first describe some intricacies in the gen- 
eral analysis of ( pTl] ) if V{x) is not finitely supported (as 
would be natural to assume, e.g., in models of cavity 
QED 1^ ) , and set the stage by giving the exact solution 
of the RWE at Q. The former (Section II A) is neces- 
sarily treated in detail, since the subsequent analysis will 
turn out to involve precisely the most subtle aspects of 
the theory. However, readers familiar with the issue, or 
only interested in results on black holes, can initially skip 



over this part. In Section [II, the RWE is studied on the 
left (near r — 1). It is found that the SUSY genera- 
tor (pT^), ( |1.7D is not outgoing into the horizon, as one 
might have been tempted to assume; the consequences 
for SUSY are discussed. A logical question — which will 
be answered in the negative — is whether any of the physi- 
cally relevant solutions encountered in the analysis would 
be distinguished by its behavior near the black-hole sin- 
gularity r = 0. Since this involves methods similar to 
the ones at r = 1 . and since some of the results are rele- 
vant to Section [v] ( there even exists a logical connection 
back to Section |l|, but that can be established only in 
Section [v| ) , the analysis near r = is presented in Sec- 
tion 1^. In Section 0, the RWE is studied on the right 
(near r = oo); combination of the results of Sections [II- 
will then prove the above claims. In Section |v3, our re- 



sults are compared and combined with the Leaver-series 
solution |20[ for the outgoing wave to infinity. 



Section VII discusses the special modes of Kerr holes; 
for a > 0, these are all found to be TTMs, with both 
a TTMl and a TTMr (depending on the sign of the 
spin) present in the lower half cj-plane for each {£, to, a). 
Hence, for the TTMr this property is not conserved in 
the Schwarzschild limit (since neither the RWE nor the 
ZE has one, as has been mentioned above), and the rea- 



2 



son for this is elucidated. As a consequence of the analy- 
sis, it is found that the Schwarzschild special frequency is 
a limit point not only for the aforementioned TTMs, but 
for a multiplet of QNMs as well. Finally, Section VIII 



contains concluding remarks, discussing both the litera- 
ture cited above and some remaining questions. 



II. PRELIMINARIES 

A. Outgoing waves for long-range potentials 

W hi le it is well known that the tails of the potentials 
(1.3), (1.4) make the determination of their modes very 
hard numerically unless at least | Im a; | < | Re a; | , not 
all authors have sufficiently emphasized several difficul- 
ties of principle which occur on the NIA (so-called zero- 
modes), as already illustrated by the surprising coinci- 
dence QNM^TTMl at the end of Section |. Define the 
out goin g waves to the left (/) and right (g) as solutions 
of (LI) obeying 



fix- 



'OO, LU ) 



1 • e" 



g(a;^oo,tj)-l-e*'^^ (2.1) 



for Imw > and further by analytic continuation; these 
are the functions figuring in the retarded Green's function 
G^{x,y;Lu). The continuation conserves the asymptotic 
forms {'2A) for those uj for which the functions are finite, 
as long as |7r/2 — arg(a;)| < 37r/2. A solution is said to be 
incoming at lo iff it is outgoing at —lo. The situation, as 
recently examined by us in ||l^ in the context of SUSY, 
now is as follows. 

First, it could happen that some modes are of higher 
order, that is, under a generic perturbation they would 
be split up into several modes rather than merely shifted 
1^,^. However, since it will not occur in the present 
work this possibility will not be pursued here. 

Second, in general the outgoing waves will have a 
branch point at w = if the tail is stronger than expo- 
nential. Eqs. ( |1.2[ )-( |l^ show that this happens for g but 
not for /; now the NIA is the only choice for the branch 
cut which respects the symmetry g{—uj*) — gico)*. On 
the cut, one defines gi(^,.-f{uj) — limej^o5(w T e), and 
Sg{Lj) = gr{cd) — gi{to). Since (2.1) shows that the in- 



creasing parts of 51/]. (if present at all, see below) cancel 
in 6g, and since for each non-real uj the KGE (1.1) has a 
unique decreasing solution, one must have proportional- 
ity 



Sg{x,u) = a{uj)g{x, -w) 



(2.2) 



with a purely imaginary p3| . For each uj' one further 
defines a discontinuity index fi by Sg{uj) ^ (lo — uj')^ . 
Generically one expects = while no cut at all would 
amount to /i = oo, but at isolated frequencies, corre- 
sponding to the zeros of a, positive integer values for fi 



are possible as well. For further developments, see Ap- 
pendix 

Third, if the potentials are not oscillating then frequen- 
are the only candidates for anomalous 



cies LO 



points, where 



f(c.) 



uj'Y 



(2.3) 



with f(cij) finite and nonzero near oj = oj' and with the 
anomalous-point index v > I (of course the same can hap- 
pen for g, but this will not occur in the present paper). 
While these divergences in the outgoing wave are reason- 
ably well known, up to now it apparently has not been 
realized that they are a mere artifact of the normaliza- 
tion of /, so that for instance in G^{uj') they would cancel 
against the corresponding divergence in the normalizing 
Wronskian Rather, the true outgoing solution to 

the left at the anomalous point is Since {2A) and 

dJ) show that f(w) - {uj-uj'Ye-''^=', it follows that the 
increasing part of f{u>') vanishes. Hence, f{uj') is actu- 
ally decreasing for x — > — oo and thus coincides with the 
incoming wave: f(w') oc /(— w'). Because of this coun- 
terintuitive behavior, f(cij') is called anomalous outgoing. 
It should be noted that lim^^-^i^' lmix~,-oo — °° 

but lima;^_oo limi^^o;' |f(a;,w)| = 0, severely hampering 
any brute-force numerical analysis. However, ^(x,lo) is 
completely regular near uj' in any finite region of x. 

Physically, the above merely means that at the anoma- 
lous point lo' , the small (typically exponential) tail scat- 
ters so strongly that the outgoing wave is completely dif- 
ferent from the free outgoing wave; that is, one has a 
resonance phenomenon in 1 -f 1 dimensions. Technically, 
the derivation just given is almost identical to the one 
of the simple proportionality 8g(x,Lo) cx g(x,—u}) in the 
preceding paragraph. However, the concept of anoma- 
lous points, while central to this paper, can initially be 
confusing si nce the definitions of outgoing and incoming 
waves near (2.1) would at first sight seem to be incom- 
patible unless LO — One can, if necessary, convince 
oneself by taking Vix) — VqC'*'^ (A > 0) and calculate at 
least k Born approximations to j[uo) to see the anoma- 
lous point at lo' — —ik\/2 {k — \ being the simplest 
example). Rescaling / i— > f as in the preceding to remove 
the anomalous-point divergence, one observes that f(w') 
and f{—Lo') agree to any desired order, as they must. 
Considering the exact solutions available in this or 
in the Poschl-Teller [|l^ case, using standard identities 
for the Bessel and hypergeometric functions involved one 
arrives at the same conclusion. 

These examples also introduce the following approach: 
if V^(A^^lnz), continued from z > 0, is analytic near 
z — Q for certain A > (the assumption of a decreasing 
potential tail implying that V{z^Q) = 0), the KGE ( p^ ) 
has a regular singular point p5[ at z — in the variable 
z = e'^^, and the index equation yields the associated 
characteristic exponents simply as ±ia;/A (from now on 
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A = 1 not to burden the notation) . The Frobenius theory 
of regular singular points now distinguishes three cases 
p6t . (a) In the generic case 2ia; ^ Z, there are two inde- 
pendent generalized power-series solutions, correspond- 
ing to the Born series for the outgoing and incoming 
waves. If 2iuj £ Z, the small z*" for Imw < 0) solu- 
tion still is a generalized power series; this is the incoming 
wave /(— w). However, the large solution is of the form 



Eoo 
.7=0 ' 



+ Cfi^j ^i^z, where {cj>i} and C have 



to be calculated case by case, (b) Typically ^ ^ 0, so that 
there is no generalized power series for the large solution; 
this is exactly the case of diverging Born series discussed 
before, and hence corresponds to an anomalous point. 
Indeed, in the example of the preceding paragraph these 
were located at uj' = —ik/2. The outgoing wave, defined 
by analytic continuation, will certainly not be logarith- 
mic in z and hence should correspond to the small solu- 
tion, as has been deduced already for anomalous points 
in general, (c) It can also happen that a 'miracle' ||ic| ] 
makes C — even though 2iuj S Z. Hence, in this dou- 
bly nongeneric situation, there is a one-parameter fam- 
ily of large generalized power-series solutions, undeter- 
mined up to a multiple of the small solution. One of 
these will correspond to the analytically continued out- 
going wave, since in the coordinate x the vanishing of 
C means that higher-order corrections to V{x) ~ Vbe^ 
have conspired in a 'miraculous' cancellation of the di- 
vergence in /(w). Thus, case (c) does not correspond to 
an anomalous point. 

While the Frobenius method in itself of course is ele- 
mentary, the above recapitulation emphasizes its aspects 
when Lo is not a fixed parameter, but a variable of the the- 
ory. Let us complement the discussion by two remarks, 
one more abstract, the other more computational and 
concrete. 

The first remark is that one can analytically continue 
the solutions along a circle around z = 0, and study 
the so-called monodromy map this generates in the two- 
dimensional solution space For generalized power 
series the effect follows from the leading coefficient as 
^±iuj ^ (^^^2^^■^±^u ^ ^±^u^T27ruJ^ Heucc, in the generic 

case (a), the monodromy map is a diagonalizable ma- 
trix with distinct eigenvalues e^^^'^. On the other hand, 
in the anomalous-point case (b) this map is a nontrivial 
Jordan block (e.g., |^), because applying it to any large 
solution yields an image which is not even proportional, 
as the continuation carries one onto a different branch of 
the logarithm. The missing eigenvector precisely corre- 
sponds to the missing preferred large solution, since both 
the incoming and outgoing waves are small. This also 
shows that, if one were to define a preferred 'interesting' 
large solution by other means, such a hypothetical notion 
of 'interesting' cannot be invariant under monodromy, 
which limits its relevance. Finally, for a 'miracle' (c) the 
monodromy map becomes either —1 (for half-integer icu) 
or (for integer iu), so that any solution is an eigen- 



vector. Hence, in the latter case, at the immediate level 
only the incoming solution is distinguished as the small 
one, while the preferred increasing solution — the outgo- 
ing wave — must be determined by other means. 

The second remark concerns precisely the calculation 
of such miraculous outgoing waves. Generically the co- 
efficients of f{uj) — J^JLo CjZ^~"^'^ follow from the recur- 
rence relation j{j — 2iuj)cj — J2p=i^p'^j-p^ with V{z) = 
E^i VpZ^ and Co = 1. If uj = uj' ^ -ik/2 [k e Z), the 
brute- force approach is to calculate Ck{uj) from this recur- 
rence for unevaluated uj, and finally investigate the limit 
UJ ^ uj' . However, even using computer algebra such op- 
erations on (big) rational functions are less transparent 
than those on (large) integers only. Of course, it is stan- 
dard to verify whether a miracle occurs by checking if 
X]p=i ypCk-p(uj') vanishes. If it does, in this way one has 
verified that a large outgoing wave exists. However, even 
to determine which value of Cfc corresponds to J{uj') brute 
force is not needed. Namely, using de I'Hospital's rule 
one finds Cfe(w') = {-2ik)-'^ Y!i,=i^p(^^c.k-p\uj' , where 
for 1 < p < k — 1 the numbers di^Cp\uj' follow from the 
w-derivative of the recursion relation, which in the given 
range of p can never lead to any divergences. Finally, the 
Cj>k+i are calculated from the recursion relation with- 
out further ado. For the RWE, in the following sections 
we will obtain complete analytical control, and a double 
check (which we have performed) using brute-force com- 
puter algebra is still feasible for £ = 2,3. However, in a 



more complicated situation in Section VII B we shall use 



an approach analogous to the above to calculate (7.24). 

When it applies (as it does for black holes), Frobe- 
nius theory thus yields complete information on anoma- 
lous points and miracles. However, other cases are 
also possible, such as V{x) ~ xe^^ (readily handled as 
d\[V{x) e^^], yielding an index i/ = 2 jl^), and the re- 
sponsible exponential tail could even be buried beneath, 
e.g., an algebraic one; see further Section VA. 



B. Exact solution of the Regge— Wheeler equation 



To avoid the transcendental ( |l.2| ) we write the RWE 
in terms of r, 



[r\r-lfdl + r{r-l)dr 
- (r-l){2(n+l)r - 3} + LU^r'^](j) = , 



while (2T) translates to f{r,u}) 
g{r,uj) ~ {re^'Y^. At = Vt?, integration of (|l.5D-(| 



(2.4) 

e(r-l)]-''^ and 



yields as the 7 = 71 = case in the exact solution |l4] 

(2.5) 



C(r,7) = [(r-l)e'l-^/2^^^ [1+7/(^)1 , 



dt 
t- 1 



{2nt + ?>Y ' 



(2.6) 
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readily found by varying the constant in ^i. The re- 
maining quadrature seems to make ^ (r, 7) singular near 
2nr + 3 = 0; however, in fact this cannot happen since 



the RWE (2.4) is regular there, so the singularity must 
cancel because n and N in (2.6) are related as in (1.6). 
Motivated by these remarks, one finds that the integral 
actually is elementary, viz., 



(2n-l)(n+l)(Ar-2) 



N-2 
3=0 



[Nil-rW 



e^''(r-l)^-i[2nr2 - (2n+3)r + 6] 
4n2iV(2nr + 3) 



(2.7) 



Various choices of the free parameter 7 now will cor- 
respond to the various physically interesting solutions. 
We have already defined ^1 with 71 = as the SUSY 
generator (L7). Since incoming waves are readily and 
uniquely characterized by being decreasing, their deter- 
mination is immediate. In particular we s ee that ^1 itself 
is incoming from infinity, and from (2.5) we get in de- 
tail ^1 — 2ng{~il). Similarly, the form (2.5) shows that 



^2 = lim^^oo ^(7)/7 = /(-^)/[(2?^+3)iV] is incoming 
from the horizon, i.e., formally 72 = 00. Further, we 
define ^3 as the outgoing wave into the horizon cx /(r2), 
^4 as the unique small solution near r — 0, and ^5 as the 
outgoing wave to infinity oc .g(r2). It may be helpful to 
note that for 1 < j < 5, £,j is evaluated in Section j. 



For general uj, the RWE (2.4) has regular singular 
points at r = 1 and r — 0, and an irregular singular 
point at r = 00 |l5|. We shall successively study these 
three singularities (in fact the behaviors of the solutions 
at each of them will turn out to be related); upon com- 
bination, the results will completely elucidate the situa- 
tion at the special frequency and in particular prove the 
claims made in the Introduction. 



III. NEAR THE HORIZON 



At the horizon r 
described in Section 



1, one has exactly the situation 



II A 



(with A = 1), so that the index 
equation has roots ztiut. Hence, for 2iuj — k ^ Z\{0}, 
typically the large (r— l)^!"^!) solution is expected 
to contain a contribution ~ (r— l)!"^' ln(r— 1) so that 
Lu — —ik/2 is an anomalous point, as is indeed read- 
ily verified for the lowest few k. However, one now sees 
at once that, for any £, lu = Q is a doubly nongeneric, 
'miraculous' frequency, for the general solution ^(r, 7) of 
(2.5) contains no log-terms at all. Hence, the Born series 
is finite and convergent, and the outgoing solution /(fi) 
(up to this stage: whatever it may turn out to be) is large 
so that Lu — Q is not an anomalous point. If it seems sur- 
prising that, e.g., the simple form of ^1 alone can lead to 
such a nontrivial conclusion, one can of course for a given 



£ evaluate the first TV — 1 terms of the series for /(fi) to 
see the pertinent cancellation (Section II A) occur. 

In contrast to the incoming wave, in the presence of po- 
tential tails the outgoing wave is no longer characterized 
by its leading asymptotic behavior, nor does the solution 
( p.5D at a single frequency allow the continuation proce- 
dure stipulated below (2.1). An elegant indirect method 



to find the outgoing wave for r | 1 is as follows: £, is out- 
going iff dc^[{r—iy^£] contains no log-terms. Namely, in 
the general superposition di^[{r—iy'^{Af{uj)+Bf{—iLj)}] 
only the second term in square brackets contains u in an 
exponent, so that these log-terms come from the incom- 
ing part of ^. Writing things out using the variation-of- 
constant solution for d^^S, |2^, the criterion becomes that 
[r/(r— l)]^^(r) (the first factor is a Jacobi an) should not 
contain an (r— l)^^-term. Substituting ( p.5|) into this 
nonlinear function of 1^, one finds that this happens iff 



7 = 73 



2 ^ 



28-11093 



-|-3(2n+3) 



for 



N 



N+l 



Nl 



-N 



(3.1) 



(3.2) 



in particular, 73 is a sum of negative terms and thus 
is nonzero for any £. Hence, defining ^3 = ^(73) = 
{2n+3)f{n), for all £ one has £,1 <^ £,3 <^ £,2, where the for- 
mer inequality means that the RWE has no TTMr at fl. 
Summarizing, we are led to the striking conclusion that 
in spite of its appealing simple form the SUSY generator 
^1 is not outgoing into the horizon. 



The special case (3.2) for the 'magic number' 73 can be 



and has been readily verified by using computer algebra 
to do an eighth-order Born approximation with uneval- 
uate d LP, or its refinement described at the end of Sec- 
tion II A . However, for the next-simplest case £ = 3 one 
already needs fortieth order Despite this computa- 
tional demand, these series methods are a reliable way of 
computing f{uj), the logarithmic derivative of which can 
be compared to the one of g(— tj), the latter being stably 
obtained by integrating the RWE down from large r. In 
this way we have verified that for £ ^ 2 there is also no 
TTMr near il on the NIA (say between the nearest two 
anomalous points, although the rest of the semi- axis is 
readily checked as well). In fact, for any £ one readily 
convinces oneself analytically that the real zero of f{x, il) 
(cf. the asymptotics of ^3 for x ±00) is stable under 
small imaginary perturbations of H.. Since for a repul- 
sive potential such a zero means that /(w ~ Q) cannot 
simultaneously by incoming from the right, the claim of 
]l7t on the presence of a TTMr is refuted. 

The above is enough to classify the RWEk->ZE SUSY- 
transform on the left. Namely, has turned out not 
to be anomalous, while the generator ^1 is mixed (i.e., a 
nontrivial superposition of incoming and outgoing). This 
means that the transformation is of category (b2) in the 
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terminology of Appendix B3 in [Q, the result quoted 
in Section VI of that paper. We then get the prediction 
that the 'miracle' will not repeat itself for the ZE, where 
on the left should be an anomalous point with index 
i> = 1, being the only solution with no log-terms 
besides, it follows from the general theory that cx f (fi) 
is anomalous outgoing. Independent of the precise situa- 
tion on the right, this means that for the ZE a TTMr at 
n would simultaneously be an NM, which cannot occur. 
Hence, also the ZE has no TTM^ at Q,. 

Let us conclude this section by pointing out two in- 
teresting perspectives on SUSY offered by the preceding. 
In the first place, in we obtained a very straightfor- 
ward proof of ^1 not being outgoing into the horizon — a 
property referring to the RWE only — by looking at the 
anomalous point of the ZE. In the second place, in Sec- 
tion II A it has been argued that this very anomalous 
point is a nuisance numerically near f2. However, in this 
section we have seen that the corresponding anomalous 
point is absent in the RWE, and any numerical result 
obtained there can subsequently be SUSY-transformed 
back to the ZE. Thus, SUSY can be used both as an an- 
alytical tool for the RWE and as a numerical tool for 
the ZE. The theory of ||l^ shows that the latter applica- 
tion actually has some generality, namely, as long as the 
anomalous-outgoing wave is nodeless (as it has to be for 
repulsive V, cf. pSj). 



IV. NEAR THE BLACK-HOLE SINGULARITY 

The index equation now gives the characteristic expo- 
nents pi — 3, p2 — —1, so that pi — P2 G Z. While 
this is still widely recognized (e.g., (20j), previous works 
do not seem to carry out the associated series expansion 
for a more precise determination of the character of the 
singularity, i.e., logarithmic/anomalous (case (b) of Sec- 
tion 



HA) or miraculous (case (c)). The above shows that 



this determination in fact is decisively simpler than for 
r = 1, since the expansion is only needed to fourth order 
for any £ and uj. By again merely looking at the form of 
(2.5) (cf. Section HI), one sees that at least ioi u; — H. a 
miracle happens at r = 0, since the general solution then 
contains no Inr-terms. Indeed, actually carrying out the 
fourth-order expansion readily shows that a power-series 
solution ~ r~^ exists if and only if uj^ = 17^. 

Let us now define ^4 as the unique solution which is 
finite near r = 0: 



^4(0 - ( — 



Ar/2 



2nr + 3 f dt [(t-l)e*]^t3 



r Jo t-1 (2ni + 3)2 
= ei(r)[-/(0)+/(r)] 

= -/(0)ei(r)+6(0 , 



(4.1) 



(4.2) 



so that 74 = — /(O) Obviously ^4 9^ ^1, since their 
r — > 00 asymptotics differ. Also ^4 9^ C2 since /(O) 7^ 0, 



all factors of the integrand having a fixed sign on (0, 1). 

Could ^4 ^3 = ^1 -I- 73^2, that is, could the un- 
expected result ^3 7^ ^1 for the outgoing wave into the 
horizon perhaps be distinguished by its behavior near 
the black-hole singularity? This is true iff /(O) — — 7;^^- 
We will now prove that this cannot happen; for the re- 
quired bound, (^) is less convenient than direct es- 
timation of the integral. Because N is even one has 
/(O) > 0, and /(O) < y^dt{l ~ t)^'^e^H^; use 

J^dt{l^t)'^e^' = (fc!/A^^-+i)[e^ - E.to^Vj!], do- 
ing the terms with N < j < N + 2 exactly. In the 
remainder use 

- E"=o ^Vj! > 0, then the con- 
tribution of this remainder is seen to be negative, so 
it can be omitted in the upper bound. After this, 
the inequality becomes /(O) < 2(A^-f 1)/9A^^ < 1/37V^ 
For 1 73 1 use YI^^q N^/j\ < e^, 2n + 3 < 3^77, and 



m > 2A^^+i/2e-^, yielding 
I73I <9N 



(4.3) 



Combination now gives /(0)|73| < 3/A^ < 1, for N > 8. 
Hence, ^1 through ^4 are all different. 



V. NEAR r = 00; THE MODE SITUATION 
A. A continuation theorem 

Let us first derive a general result on the KGE (|1 . l|) : 

if V{x) is analytic for x > Xq (say), and if its continua- 
tion from those x has V(x) = 0{x^^^^) for \ argx\ < (3 
and some /3,7y > 0, then glx,uj) - 1 • e^'^^ if (A) 
\argx\ < 13, (B) \tt/2 - argw| < 7r/2 -I- (3, and (C) 
\t: /2 — avgujx\ < 37r/2. In particular, g{uj) has no anoma- 
lous points for \t: 12 — avgu}\ < min(7r/2 -f /3, 37r/2). (The 
counterpart for / is obvious but not needed here.) 



Namely, for oj obeying (B) the continuation of (2.1) 
can be first carried out in parts of the x-plane where it 
merely amounts to selecting and normalizing the unique 
small solution, which can never lead to any divergence. 
This proves the result for |7r/2 — argwxj < 7r/2, which in 
a second step can be extended at fixed u by trivially solv- 
ing (1.1) in a region where V is negligible; condition (C) 
(implied by (A) and (B) if /3 < 7r/2) arises because the 
asymptotic form for g cannot be continued in this way 
from a region where it dominates across an anti-Stokes 
line. 

In applications, the second step would typically take 
the form of continuation back to the physical region 
X > xq. The rationale for the two-step procedure is seen 
to be that continuations in x are easier than those in w, 
since for the former one has the linear ordinary differ- 



ential equation (1.1) at one's disposal. For an example, 
calculating a few terms in the Born series shows the the- 
orem in action for V{x>xq) — "^^^ 
has f3 = arctan(A/|/i|). 



e cos{px) fUl, which 
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B. Application to the Regge— Wheeler equation 

Turning to the RWE, the above procedure is best car- 
ried out in terms of r. Since V{r) vanishes if |r| — > oo, one 
has /3 = oo so that g has no anomalous points, its asymp- 
totics being as below ( ^.4| ) for \it/2 — argwrj < Sir/ 2 
[^,|^. This should now be compared with the asymp- 
totics of ^, for which (2.5) and ( |2.7D show that 



^(r,7) - 2n 



1+7 



(2n-l)(n+l)(A^-2)!e^ 



7 

H [re 

2nN 



(5.1) 



in the entire r-plane. In particular, for both g\ (which 
has argfi = 37r/2) and gj. (with argil = — 7r/2) one can 
make the comparison with (5J_) in the region Rer < 0, 
where the asymptotics are not dominated by the second 
term on the r.h.s., as well as for the original a; > 0. This 
of course reproduces that ^(7) is incoming from infinity 
for 7 = 0, but also shows that it is outgoing for both 
branches of g{uj) iff 



iriN^+^e-^ 



1 ^15 



{I - 2n){n + 1){N - 2)\ 



(5.2) 



in which case ^5 = ^(75) = (75/2nA^)(7i/r(il), in accor- 
dance with the general statement (see iQ, Appendix B2) 
that if g\ and g^ coincide they must be real |35[ . Hence, 
the discontinuity index /i > 1 at fi. 

This statement on /i (dealing with analytic continua- 
tions in w and with r — s- 00) has thus been related to 
the 'miracles' near r = 1 and r = (dealing with con- 
tinuations for small r at a fixed w), since the latter make 
that C(r, 7) has the same asymptotics for argr < — 7r/2 
and argr > it/2. Note that for a qualitatively similar but 
otherwise arbitrary potential, the miracle is a 'leftist' and 
the vanishing of a{Lo) a 'rightist' concept, so that the two 
would be completely unrelated. Of course, the RW po- 
tential is anything but arbitrary , being of the form 
V{x ) = W'^{x) — d^W{x) + Vi^, as follows at once from 



(1.2) , (1.3), and (L5) or indeed from SUSY in general. 
Obviously ^1 9^ 9^ ^2: where the former inequality 

means that Vl indeed is not anomalous for the RWE on 
the right, and the latter that the RWE has no TTMi^ 
at Vl. Two questions remain, namely, whether possibly 
75 = 73 or 75 = 74; the answer is negative both times. 
For the latter (Section |^ showed that the outgoing wave 
into the horizon was not distinguished near r — Q, but 
could perhaps the outgoing wave to infinity be?), this 
is immediate since (2.7) for /(O) — —74"^ shows that 
7^^ — ^'^^ is a sum of positive terms, hence nonzero. 
For the former, using methods as above ( [4. 3D (in partic- 
ular, iV! < V^N^+^l^e-^) leads to I75I > (Vih)N^. 
If £ > 3 =^ iV > 40, this imphes I75I > 97V > I73I (cf. 

(4.3) ), while for the remaining case i = 2 the inequality 



?5 9^ 3.nd the RWE has no QNM at fi, even though the 
positive discontinuity index /x > 1 in itself would have al- 
lowed this. Thus, the five solutions we have identified — 
while of course not independent — are all different, and 
the surprise is rather that there are not six of them, i.e., 
that gi{il) = gri^), given that the outgoing wave for the 
Schwarzschild black hole definitely has a branch cut in 
general (see and Section [VD|). 



C. SUSY and the Zerilli equation 



The mere absence of anomalous points in g for the 
RWE (cf. above ^)) already implies that the SUSY 
generated by ^1 is of type (al) on the right in the ter- 
minology of The general theory then shows that 
^1 oc gi{Q) — gr{^)- Since we have shown in Section III 
that it is also anomalous outgoing on the left, ^1 thus is a 
QNM^TTMi^ of the ZE (of course, this once more shows 
that the ZE has no TTMr at 9). 

Because /i > 1 already in the original system, the 
equality of gxjj. in the above can immediately be sharp- 
ened to /i > 2 by the elementary SUSY relationship 



a(ijj) 



■ LO 



a(uj) 



(5.3) 



However, to determine the orders of the modes one needs 
the fuU apparatus of [0. Its Eqs. (B6), (B7) show that 
the QNM Wronskian Jq(w) is regular near fi, while the 
TTM Wronskian Jt(ijj) has a first-order zero there and 
a simple pole at —VL. Accounting for the anomalous 
point of /, this means that hath the QNM and the TTMi^ 
of the ZE furnished by ^1 are simple — intuitively plau- 
sible, since the RWE has no QNMs/TTMls to double 
the ones generated by SUSY. The proportionalities in- 
volved are elucidated as ^1 = ^2 =^3/73 = ^4 ='^5/75 = 
f{-Vi)/{2n+2,) = -i{2n+S)-f3^1{n) = g{n)/{2n). 

Incidentally, since ydrj— +cxd) vanishes just as well as 
the RW potential, g has the same asymptotics as estab- 
lished above (5.1) for g. Analogously continuing ^i(r) to 



175 1 



e-**2^7l35 > I73I (cf. (3.2)) is immediate. Hence, 



Re r < then reproduces that this function is outgoing 
on the right, without having to calculate the general solu- 
tion of the ZE at fl^ and independent of SUSY. This 
also shows why the ZE having a{fl) — in the absence 
of a 'miracle' does not contradict the discussion below 
( ^.2|) , relating the two. Namely, in this case g, being a 
TTMl, is small and hence single- valued near the branch 
point of the general solution. 

For readers mainly interested in results, it should be 
pointed out that this subsection inter alia proves that 
the ZE retarded Green's function has a pole at w = 17, 
the statement of which is completely independent of the 
intricacies of anomalous outgoing waves. 
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D. Discontinuity indices 

The above has estabhshed the mode situation at the 
special frequency. It will now be shown that the derived 
inequalities for the discontinuity indices are in fact real- 
ized as equalities, i.e., that one has fj, = 1, which implies 
/i = 2 by (5^). To this end, define gi{uj) = d^g{uo) 
pl| and make the variation-of-constant Ansatz gi{^) = 
g[n)h[n) [0. One finds 

r t 

d^h{x,^) ^ ~2Vlg-^{r,n) i dt g^it,n) , (5.4) 



t-1- 

and to get the correct asymptotics dxh{x, f2) — + i also 
if r is continued to the left half-plane, one should set 
To = —CO. However, for gii{r,n) {gii{r,il)) the asymp- 
totics are imposed if r is continued into the second (third) 
quadrant, so that for r > one must put t t + ie 
{t 1-^ t — ie) in the integrand of ( |3.4[ ). Evaluating the 
residues at t = 1 and t = 0. one finds in detail 



t-i 



^(73-75) + 9( — 



1 



(see Section II A for the definition of S). In particu- 



that 



lar, since it has already been shown in Section 
1 75 1 > 1 73 1, the first term in brackets of (5.5) (which 
is due to the pole at i = 1, and hence anticipated to 
be oc 73 — 75) is positive as well as the second so that 
Sgi{n) ^ and = 1 11]. Eqs. (U) and dJ) wiU have 
an interesting application in Section VIII. 



VI. LEAVER SERIES 

In we already touched upon the relation of this 
work to the first Leaver series 



9 ^ \IUJ OC 



3=0 



r- 1 



(6.1) 



However, at least as interesting a connection exists with 
the second Leaver series. 



where M4 is the series to four terms of the regular con- 
fluent hypergeometric function M, and ip = F'/F. 

To see the connection, let 2iuj — L = 1,2,..., so 
that the Pochhammer symbol reduces to (1 — L)j = 
{-y{L-l)l/{L-j-l)l for < j < L - 1, and zero 
for j > L. Typically this zero is compensated by the 
anomalous-point divergence in Cj, so that the series still 
involves C/(5,5, z) = e^F(— 4, z) and hig her, which have 
logarithmic branch points at z = (cf. ( |6.3| )). However, 
if = LUai is a 'miraculous' frequency, the sum in ( |6.2| ) 
truncates: 



7(r, Win) oc 



r- 1 



L/2 L-1 



X U{3+j-L,5,-Lr) 



(L-J-l) 



(6.4) 



in which only those Cj's enter which are the same 
for all solutions which are large near r = 1, so that 
in particular for L = N one has Cj = (1/j!) x 



rfiL=oexp{A^/(u-l)}[2n(l-u)^ + 3(l-u)^+i] (comb: 



( p.5[ ) and (|6T|)). The functions U involved in 
evaluated from (|6.3D as Q 



me 
are 



C/(2,5,z) = z-^ + 4z-^ + 6z' 



C/(l,5,z) = z-^ + 3z~^ 
k 



6z' 



6z 



Ui-k,5,z)^J2 



p=0 



P 



{p + b)k^p{~f+^Z^ 



(fc = 0,l,...) 



(6.5) 



Since f7(2, 5,z) and [/(l,5,z) are linearly independent 
and since cl-i and cl-2 cannot both vanish for a nonzero 
sequence {cj} because of the recursion 



- [2j2 + (2-4L)j + 2L{L-1) + 2(n+l) 
+ if - 2Lj + L2 - 4)c,_i - , 



3]cj 



(6.6) 



g{r, uj) oc r' 



X U{3+j-2iLu,5,-2iujr) 



in terms of the same cj as in (S.l), and one of the more 
'irregular' cases (integer second argument) of the irregu- 
lar confluent hypergeometric function [p9[. 



U (k, 5, z) 



1 



24F(k-4) 



M(k, 5,z)lnz 



00 / \ p 



p=0 



(5), p\ 



+ 



F(k)z' 



■M4(k-4, -3, z) , 



(6.3) 



one concludes that g(r, Wm) is a large single- valued so- 
lution near r = 0. In Section IV it has already been 

n^, so 



(6.2) established that this can only happen for ui^ 



that n is the only 'miraculous' frequency of the RWE. It 
should be emphasized that this is a nontrivial statement 
on the singularity near r = 1, and that the above proof 
involved an analysis both near r = and r = 00. 

Thus, L — N \s the only case which actually oc- 
curs, and hence the r.h.s. of (6.4) should evaluate to 
2{n+\)N-^£,r^{r) as given by (|2.5|), jO]) and (|5^), where 
the constant of proportionality (for the Cj below (3.4)) 
has been determined from the leading asymptotics for 
r — > 00. Verification amounts to comparing two polyno- 
mials of degree iV-|- 1 with rational coefficients, which has 
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been done explicitly for i = 2. For general £, the stated 



outcome of the sum in (6.4) can be proved by downward 
induction in the power of r, using the recursion (6.6) up- 
ward; however, we omit the laborious details. 

The expression (6^), ( |6.3| ) is also of interest for non- 
miraculous uj — Namely, it shows that in the 
Leaver series for bq one has 



(5[/(3+j-2it^, 5, -2itjr) 



iTrAf (3-|-j— 2?w, 5, —2iLor) 



12r(j-2itj-l) 



(6.7) 



The relevance of this latter relation is not so much that 
the regular M is easier to handle than the U function 
(which it is), but rather that it points to a computa- 
tional scheme for 5g directly and not as an exponentially 
small difference between the gi/,.. In particular, the pro- 



portionality derived in Section [I A means that one can 
evaluate a{Lo) at any convenient r > 1, and a method of 
calculation which avoids substracting nearly equal large 
numbers would enable one to look for, e.g., possible sim- 
plifications as r ^ oo. However, the properties of the 



(non-truncating) sum (6^) itself on the NIA are still un- 
der investigation [^. In its turn, this last remark also 
motivates the explicit verification of (^^) in the previous 
paragraph. 



VII. KERR HOLES 

A. Kerr-hole special modes 

Both the exact solution for the stationary metric and 
the separated wave equations for its linearized perturba- 
tions can be generalized from the spherically symmetric 
Schwarzschild hole to the axisymmetric, rotating Kerr 
hole. Presently we investigate the issues addressed in 
the preceding sections in this wider context, and start 
by establishing some notation. Expanding the pertinent 
effective scalar field as 

cx; £ 

^(u;,r,0,^)=^ ^ e""^5(0,c.)i?(r,u;) (7.1) 
in Boyer-Lindquist coordinates, the radial equation reads 



A'dr 



drP 



K"^ + is{l-2r)K 



+ iisujr — A 



^ , , sK l-2r\/ , .sK , ^ 

2A'^ — A~ / l *2A ' '^^^ ^ 



where we introduced the auxiliary quantities 



P = 

P = 0, 
(7.2) 



A = - r + 
= (r-r+)(r 

1 ±6 
r± = -IT- ' 



b= Vl - 4a2 , 
K ^ {r^ +0^)10 



(7.3) 
(7.4) 

(7.5) 
(7.6) 



^ is the 



and where P ^ A^+'/^i?. In ([tJ), < a < 
black-hole rotation, and s = ±2 is the spin of the field 
The radial field P depends on (r, oj, £, m, a) and on the 
sign of s, but this will only be indicated where necessary. 
The angular separation constant A is an eigenvalue of the 
equation for S{9, cu), cf. ( 7.22 ) below, and different eigen- 
values correspond to different angular-momentum sec- 
tors. Since A obeys the symmetries A*„(ci;) = A_„(— w*), 
%m{^) ~ A_nt(— g;), and As = A_<; (besides, A(a-', a) = 



A(cja)), from (|7J) it follows that P^{uj) = P_^{-uj*) 
and Pmsi^^) = P-m,-si—^)', an advantage of using P 
instead of R is that this makes the latter relation man- 
ifest. Moreover, there is the considerably deeper sym- 
metry between P^ and P-s furnished by the Teukolsky- 
Starobinsky identities (which follow by combining the 
finite-a generalization of ( [7.11[ ) below with its inverse for 
opposite s, see also Q). Since the special modes ( |7.8| ) 
below are in the kernel of the corresponding symmetry 
operator (cf. (7.14)) we do not give the formulas; how- 
ever, the existence of this symmetry will be used at the 



end of Section VII C 



Where the mere existence of the ordinary differential 
equation (^2) is already remarkable, this holds a for- 
tiori for its exact solvability if uj = f2(a) is a zero of the 
so-called Starobinsky constant 

S = a2(A+2)2 + 8A[5A(amw-a2w2) ^ Q{amuj+a'^uj'^)] 
+ 36cj2 + lUiamuj-a^uj^f . (7.7) 

Namely, in that case the special modes P = V are given 

by 



V 



isilr^ + Ar^ 



sQ 
i—-r 



2a^ 



|A(a^ 



A-l-2 







J 




— r_ 





24r22 

Q A(A+2) + \2[amQ.-a^V?) + ZisVt , 
with the tortoise coordinate 



X = r ^ — — ln(r— r+) 





■ ln(r— r_) 



(7.8) 
(7.9) 

(7.10) 



mapping r_|_ < r < cxi to a; e R. In view of the aforemen- 
tioned symmetries, one can take Im f2 < for consistency 
with the preceding sections. 

When assessing the physical significance of the solu- 
tions ( |7.§| ), in comparison with the RWE/ZE case one en- 
counters two technical subtleties, which are readily dealt 



with. Writing (7.2) as a KGE in the x-coordinate for 
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\/r2~-ra^P/A, one finds that the ensuing potential is (i) 
complex and frequency dependent, and (ii) tending to a 
nonzero constant if x ^ — oo, and having an 0{x~^) tail 
for X ^ oo. (i) If V = V{x,uj), besides the outgoing 
waves one also has to define incoming waves F, G 



(temporary notation), being solutions of (1.1) at uo with 
G{x^oo,ijj) ^ for Imo; < and continued analyti- 

cally from there, and analogously for F . In the frequency- 
independent case these become simply G{uj) — g{~Lo) 
etc., but in general G and g are unrelated. Fortunately, 
for ( [7. 2D one has Gms(w) = (jf_m,-s(— w). (ii) A poten- 
tial having a nonzero spatial limit can be handled by a 
shift in frequency; in particular, for (|7.2| ) it is customary 
to define Cj = oj — ma/r+. The 0{x~^) potential tails 
give rise to a power-law prefactor in the wave-function 
asymptotics, not unlike the one occurring for the RWE 
when g is considered as a function of r, cf. below (2.4); 
these power laws will cause no problems here. 

The character of the special modes on the right is now 
easily determined. With our convention on fi, it follows 
from (7^) that V2 = 'Ps=2 is decaying and hence in- 
coming from X = -|-cx), while 7'_2 is outgoing by the 
continuation argument of Section VA. Near the hori- 
zon r = r+, the characteristic exponents of (7.2) read 



Pi(2) — 1 ± (1 -I- isrj^uj /2h). For the special modes one 
thus has pi — P2 = 2 -|- isr+il/fe, with f2 a zero of (7.7). 
For a I this becomes pi — p2 ^ 2 + sN/2 e Z, but 
one does not expect the difference to keep this value for 



nonzero a. Indeed, the expansion (7.26) below shows that 
Pi — Pi is non-integer for all sufficiently small nonzero a. 
While we have not attempted a proof, it seems very un- 
likely that, in some (£, m)-sector, the complex function 
P\[a) — /02(a) would assume a real integer value exactly 
on < a < i. Thus, for nonzero rotation, the alge- 
braically special frequency n{a) is actually generic in the 
terminology of Section II A . Since ( |7.8| ) , ( 7.1C ) show that 
the modes V are generalized power series corresponding 
to p2 (equivalently, they are eigenfunctions of the mon- 
odromy map near r+ with eigenvalue e^'^*''^), this means 
that V2 is outgoing into the horizon and hence a TTM^ 
for a > 0. Similarly, V-2 is incoming from the horizon 
and thus a TTMj^ ; this follows already from its being de- 
creasing on the left, so that for 1^-2 the condition a > 
is not needed. 

For a = 0, in which case the above analysis is incon- 



clusive for 7^2, the radial equation (7.2) does not reduce 
to the RWE or ZE, but to the Bardeen-Press equation 
instead; the explicit mapping between them reads B 



P = r-^ 



iV\ 1 

W (2nr + 2± 1) isuj 

2 I r 



X <^ d , 



ISOJ 



<P (7.11) 

= A±0, (7.12) 
where the upper (lower) signs should be chosen if 0(r, uj) 



is a solution of the RWE (ZE). Calculating the Born se- 
ries for / for (Q with £ = 2 (so that A = 4, cf. ( |7^ ) 
below), at = — 4i one finds a divergence (like for the 
ZE, but here in sixth order) if s = —2 but not if s = +2 
(like for the RWE, but here the relevant cancellation oc- 
curs in tenth order), so that is anomalous (miraculous) 
for s — —2 (-|-2). However, 7^2 is a nontrivial superpo- 
sition of outgoing and incoming waves — again involving 
the 'magic number' 11093, cf. (3.2) — in marked contrast 



to the situation for any nonzero a. Thus, the unexpected 



findings of Section [II are not an artifact of the mapping 
(7.11) — a kind of w-dependent SUSY transform — but in- 



trinsic to the Schwarzschild limit, which apparently has 
to be highly s ingu l ar. To investigate this further, one 



should study ([7.2D - (|7.8| ) for small a; this is done in the 
next subsection. 



B. Schwarzschild limit 

It should be explained why the s — 2 special mode be- 
haves differently for a = than for all other a. This is 
not a consequence of a singularity in ( [7.^ ) itself, which 
in the limit tends smoothly to 



1 



-2 {r) = — \^n^ in+iy-^ + 12n^r^ + ISnr + 9] 



X e 



-Nx/2 



(7.13a) 

-p^^lir) ^ ^'^[2{n+iy - 3r2]e^-/2 , (7.13b) 



where one should note the similarity with ( |l.4D and (1.3) 
respectively, and where the quantitative relation to the 
special solutions of the RWE/ZE reads 



s2 



N 



AM - -',-V^!'l5s,-2 . 
2n 



(7.14a) 
(7.14b) 



Since it has been shown in Section VII A that is 



outgoing on the left and that ri(a) is generic for nonzero 
a, this implies 



2np(°^ 

lim lim f2(uj,a) — 

aiO^^O(a)''^ ' ' (2n+3)(iV+l) 



(7.15) 



{2n+3)N{N+l) ' 
where /2 = fs=2 has been normalized as 

/2(x^-oo,c^,a)~l•e-^'^^ (7.16) 
and where the constant of proportionality has been ob- 



tained by combining this last relation with (7.13a) and 
(7.14a). On the other hand, by the very definition of an 
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outgoing wave, taking the limits in the opposite order 
yields an outgoing result, i.e.. 



lim lim/2(cL', a) — f2(Xl, a~0) 

uj — a|.0 



(7.17) 



(2n+3)N{N+l) ' 
where the constant of proportionality is the same as in 



(7.15) because and ^3 have equal increasing parts, cf. 
To see how these noncommuting limits arise, set 



-/b 



(i.e.. 



Co = 1), upon which the Cj follow form the recursion 
ajCj + PjCj-i + 7jCj_2 + + ^jCj-i = , (7.18) 



"J - 

1] = 
5. = 



b{3-l 

+ 4r+t2)(r+Cij— i) + 6(4ir+Cij— A) 
+ 2r+cij(2r+cij+a)) + Gibcu — A 



iCj—^ ) ( 2? — 7 — 2i(2>— ^ 
b J\ b 



(7.19) 



It is now readily seen that one has noncommuting limits, 
because 

ajv+2 = iN+2)[2ima - Na^ - 2{iLj-N/2)] + ©(ma^) 
+ 0{a^) + O[{L0+iN/2)a] + 0[{uj+iN/2)^] 

(7.20) 

implies that cn+2 (and hence higher Cj) are sensitive to 
the ratio of cj + iN/2 and a. To calculate which are the 
limiting functions, one has to iterate (7.18) up to c^r+i, 
up to the error terms indicated in (7.20). This in fact 
is a calculation of the type outlined at the end of Sec- 
tion II A, except that here one has three numerical co- 



efficients in the recursion relation — the lowest-order one 
plus the leading uj + iN/2 (linear) and a (linear for m ^ 0, 
quadratic for m = 0) corrections. In terms of [E3 



Hi) 



(^2_^2)(^2_4)2 



(7.21) 



2(£_i)^3(^+i) 

the required expansion of the separation constant reads 
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A = ^(^+l)-2- 2 



+ [h(£+l)~h{e)]a^uj^ + 0[{aojf] 
+ 0[{auj)'^] for ^ = 2 , 



2\ 2 2 

zTJi a ui 



(7.22) 
(7.23) 



and for ^ = 2 some computer algebra yields 

3633559(4-iw) + 37411744ima -t- 420400640^ 



cio 



1786050(icj - 4 - ima + 4a2) 



(7.24) 



where both the numerator and the denominator are ac- 



curate up to the error terms indicated in (7.20). From 
(7.24) it is immediate that 



lim limcio = ^ 

ui^n alO 



3633559 
1786050 



„(3) 
-10 



(7.25) 



To also find the limit corresponding to ( 7.15| ), substitute 
( |7.22[ ) into ( [7. 7] ), yielding its roots as 



n{a) = 



I- ^200^ 

O(a^) 



2„2 



0(a 



i§i(65m2 
4 



18)a^ 
for £ = 2 



no - 
n.20 — 

^21 = 



CIq + Qima 
+ 0{ma^) H 
-4i -I- ^ma 
+ 0(ma^) 

3 (2^-l)(2£+3) ' 
.4n3[21£2(^+l)2 + 2£{£+l) + 12] 
* 27(2£-l)(2£-f3) 



(7.26) 



(7.27) 



which in its own right is perhaps a new result for the 
asymptotics of the special frequencies, agreeing well with 
the numerical values in |l^ for ^ = 2 and the smallest a. 
Note that ri2o vanishes for £ = 3, and has opposite signs 



for ^ = 2 and ^ > 4. Using ( |7.27D in ( |7.24| ) now yields 

69632 



lim 

alO L 



lim 



Clo = 



893025 



(7.28) 



where in particular one obtains the same limit both for 
TO ^ and TO = even though these two cases are 
quite different for (7.24) and ( [7.271 ) individually. Since 
the normalized function lim/2(ci^,a) is a solution of the 
Bardeen-Press equation, which is completely specified by 
its coefficient C7V+2, and since for ^ = 2 the values (7.25), 



(7.28) are readily checked to correspond to (7.17) and 



( 7.15 ) respectively, the above clarifies the mechanism for 
the singularity of the Schwarzschild limit, being the zero 
in aN+2- 



C. Special QNMs 

We have seen in the preceding subsection that /2(w, a) 
is wildly varying near (uj=—iN/2, a=0). Thus, if near 
this point /2 can tend both to ^+^1 and (up to 

normalization), can it also tend to After all, the 

latter is merely a linear superposition of the former two. 
Thus, one should look for solutions of 
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/ \ (5) _ /i 75 \ (1) , 75 

cio[uj,a) = c\q' = ( 1 Cio + — 

V 73 / 



73 



2166784 
' 893025 
(7.29) 



in the region of validity of (7.24). Indeed one finds that 



hniaxo limaj^L^,(a) cio(w, a) = for 

33078176 3492608 . 



Jci 4^ 1 1 tUi 

^ 700009 



41177 



(7.30) 



This result has a significant interpretation: to leading or- 
der, for uj = ujq{a) the outgoing function into the horizon 



/2 is also outgoing to infinity, so that (7.3C) represents a 
family of QNMs of \7-^ ) for {i=2, a^O), branching from 

the special frequency fl in the Schwarzschild limit. 

Clearly, one is interested in the generalization of ( 7.30| ) 
to arb itrar y £. This can be performed in closed form, 
since ( 7.24 ), obtained in the above using computer alge- 
bra, can be reconstr ucted from its known limits together 
with the expansion ( 7.2C ). In fact, considerable cancel- 
lation occurs in the answer: while for reference we give 



the generalizations of (7.28) and (7.25) as 



-■N+2 



27{2n+3){N+l){N+2) 



18n-45) /N 



N 



[from ( |7.13a| ) and ( |7.15| )] and 



.(1) 

-■N+2 



6736 



N 



{2n+3)^N{N+lf{N+2) 



(7.31) 



(7.32) 



[from ( |7.17| ) with = A+£,i + 73^+^2 and 

= 6e^/2(r-l)^/2+2/[(27i+3)(7V+l)(7V+2)] + 
0[(r- 1)^/2+3] by the explicit form (|7tlT|) for s = 2 and 
the three leading coefficients of ^2 as yielded by ( ^.5[ ), 
( ^.6|) ] respectively, also without these explicit values one 
can derive the comparatively simple final result 



75 ~ 73 



Nj5 



2ifl20l3 . 2 

ta 



2(75 - 73) 



(7.33) 



In particular, while the 0(a)-coefficient thus keeps its 
negative sign for £ > 2, the quadratic coefficient, govern- 
ing the beh avior of the special QNMs with m = and 
positive in ( 7.3C ), is seen to become negative if € > 4. 

For a powerful and instructive check on the preceding, 
we shall presently rederive the m 7^ case (i.e., only the 
first two terms on the r.h.s.) of (7.33) from the radial 
equation (7.2) for s — —2. Namely, by the Teukolsky- 
Starobinsky identities this latter equation is physically 
equivalent to (in particular, has the same QNM spec- 
trum as) the s — 2 equation as long as 7^ f^(a), a 



condition which (7.33) fulfills for a ^ 0. To do so, we in- 
vestigate in which direction in the (w, a)-space the QNM 



( 7.13b| ) can be perturbed so that it continues to satisfy 
outgoing boundary conditions. In its turn, this calcula- 
tion buttresses our interpretation of V^l (and |i for the 
ZE) as bona fide QNMs even though they are decreasing 
on the left; however, this last property will be technically 
important in the calculation below. 

Thus, let Pq -2{a) be the branch of special QNMs, with 
u; = LUq{a) ~ ~iN/2 + ujqima + 0{a'^) and Pq^_2(0) — 

T'l^^jand set (iaPq_„2(a)|a=o = imV'^^^hq. Differentiation 
of O) yields 



drhq{r) 



pif (r) J 



-Bit) , (7.34) 



B{t) = NUq^t^ + ALOq^t^ + ( y-6cOql ]t^ 



4( ^-n-— It -2 



(7.35) 



Our limitation to linear perturbations in a enables one to 
set 6 = 1 + 0(0^) etc. in the derivation of (7.34) and 



(7.35), without which these expressions would have been 
considerably more complicated. For Pq,_2 to be outgo- 
ing to infinity, on has to set J dt = J_^dt in (7.34), by 



a repetition of the arguments in Section Since the 
integrand is entire, one does not need to specify a path; 
physically this means that our lowest-order calculation is 
free from branch-cut complications. On the other hand, 
for outgoing waves into the horizon, one sees from the 
Born series that drh{r) should have a residue 1 — Wqi at 
r — 1. Writing f^dt ~ Kq + J^dt and expanding the in- 
tegrand near t — 1, the residue criterion is seen to imply 



9c 



(1 



(7.36) 



with c being the residue at r = 1 of 



-Nr 



2ne~^'^[{N-n)r+n+l] 



(r-l)^-ir3[2(n+l)r - 3]2 9(r-l)^+2r 
e-^''[2n{'in+iy + 3(l-2n)r - 3] 



— dr 



18(r-l)^+ir2[2(n-hl)r-3] 



(7.37) 



that is, c — —73/27. For the perturbation to be outgoing 
to both sides, i.e., to represent a QNM, we now have to 
solve ujqi from 



1 -p^oru) 



(7.38) 



where a tedious but elementary calculation yields the 
r.h.s. as 8n'^{iu)qi—2nN)/^]^. Combination with (7.3(;) 



now gives Wqi = (75- 3«^^73)/(75-73), i-e., consistent 
with (7.32), which is what we set out to show. 
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VIII. DISCUSSION 



While over the years the singular nature of the extreme 
Kerr limit a 1 for rotating black holes has been increas- 
ingly well understood, the present work shows that also 
the Schwarzschild limit a | is highly delicate. Besides 
thus pointing out the problems, considerable progress has 
also been made towards their solution, culminating in 
( [7.33 ) for the special QNMs, involving a perhaps beauti- 
ful combination of 73 (Section [II), 75 (Section and 
^^ii^^20 (Section VIIB). However, besides this techni- 
cal control, it would be desirable to also have a physi- 
cal understanding. In particular, the total-transmission 
property of the 'algebraically special' solutions is related 
in to the vanishing of the Weyl scalars or ^'4 as 
these occur in the Newman-Penrose formalism of gen- 
eral relativity, an aspect outside this article's scope. At 
present it is not yet clear why the argument fails for the 
Schwarzschild hole. 

Thus, a Newman-Penrose analysis would be a welcome 
complement to our wave-equation study on the concep- 
tual side. On the practical side — the calculation of black- 
hole spectra — the prospects for numerical progress on or 
near the NIA |^l| seem good, even though the present 
work refutes essentially all numerical claims about the 
nature of the RWE at Q made previously. First, we have 
proved that the QNM fi' reported in cannot possi- 
bly coincide with fl. A spurious fi' could for instance be 
found if the numerical code ad hoc comp ensa tes for an 
anomalous-point divergence in the Cj of ( |6.l|) which at 

actually is abs ent; however, comparison of the results 
of Section VII C with the data in in fact suggests 



the existence of ft' , see below. If an imaginary fl' would 
be distinct from, but numerically very close to (even 
though it is doubtful whether the data of |15| itself sup- 
port this conclusion, cf. the Introduction), this would 
imply that a{Lu) has two very nearby zeros, which can be 
investigated numerically. 

Second, the analysis of |]l7| has some serious problems. 
While the paper does point out (below Table 1) that the 
pattern of (anti-)Stokcs lines near the horizon central to 
its analysis actually collapses on the imaginary w-axis — 
where its data are reported — the matter is not pursued 
to conclude that on that axis analytic continuations in r 
are not suitable to find outgoing waves to start with . 
Rather, on that axis the generic/anomalous/miraculous 
distinction is crucial, but the latter cannot be treated 
by the WKB ideas used in ||l^, which are more appro- 
priate for studying |a;| ^ 00 than for numerics at a fi- 
nite uj. Indeed, Refs. find that (for Rclu ^ 0) 
continuations in r can be used as a numerical tool inde- 
pendently of any WKB-type approximations. Remark- 
ably, in each angular-momentum sector (£, m), of the two 
one-parameter families — labeled by (a, s=±2) — of TTMs 
claimed in (l4|, the only one which has been comfirmed 



numerically (the TTMr at a = 0) thus is the one which 
in fact does not exist! 

Third, it has been found in Section VII C| that in each 
^-sector the Schwarzschild special frequency il{0) splits 
into one m-multiplet of special QNMs as the rotation 
is increased from zero. In particular, only for to < 
do we predict these QNMs to branch into the fourth 
cj-quadrant, in contrast with the claim in [p^ . How- 
ever, it should be noted that the coefficients in ( 7.30| ) 
are quite large, so that the Schwarzschild limit is ap- 
proached only for extremely small a; as yet no numerical 
data exist for this regime. In fact, it has been shown in 
Section VH that for the s — —2 radial equation, in this 



limit the special QNM and TTMl become identical (like 
for the ZE), while for s = 2 the single (in each {£,m)- 
sector) QNM and TTMr cancel each other (like for the 
RWE). Thus, while ||l8| correctly cautions for subtleties 
of the Schwarzschild limit, its precise nature cannot be 
elucidated from numerics alone. The preceding does not 
account for the ninth branches of (£=2,to=1,2) Kerr 
QNMs which reports for moderately small a, and 
which have slightly smaller jlmwl than their apparent 
TO — —1,-2 counterparts, the latter being accounted for 
by ( 7.3C| ) . Since in the Schwarzschild limit all Kerr QNMs 
have to merge as TO-multiplets by spherical symmetry, 
this strongly suggests that the former QNMs merge with 
their to < mirror images into an imaginary O', with 
^ 1^^(0)1 |4^. The TO = branches emerging from 
r2(0) and fi' then supposedly move towards each other 
along the NIA, until they collide and subsequently move 
away from this axis. This scenario is consistent with 
t he si gn of the 0(a^)-term in ( 7.30 ); the remark below 
( [r33|) then hints that \n'\ > \n{0)\ for i > 4. Addi- 



tional numerical work is needed to explore these fasci- 
nating possibilities; also, for a general understanding of 
the Schwarzschild limit it is unfortunate that the figures 
in do not show some Kerr modes branching from the 
high-damping series of Schwarzschild QNMs. 

A detailed computation of a(uj) as proposed in relation 
to jl^ would also be of more general interest. Namely, in 
Section we have proved that ft is the only 'miraculous' 
frequency of the RWE, and the question suggests itself 
whether it likewise is the only zero of the discontinuity 
function a. It turns out that this cannot be the case in 
general. For small lo, one can use the Born analysis of 
[Q. True to form, the Schwarzschild potentials are an 
exceptional case (a major conclusion of |Q) in that the 
leading-order contribution to a(w) expected generically 
vanishes. Fortunately the actual leading-order behavior 
still follows from the first Born approximation, as 



a{uj) ^ (-)^+i27rw + h.o.t. 



(8.1) 



incidentally for the ZE as well as for the RWE as required 



by (5.3). Eq. (B.l) can now be compared to the upshot 
of Section VD, viz.. 
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d(iuj) 

> . 



27V 



73 - 75 
7| 



,74 75 



•2) 



In both calculations one finds that the x-dependence in 
a{uj) = 5g{x,ijj)/ g{x,—uj) cancels, as it must. Compari- 
son of (|^ and (lil) now shows that the parity of zeros 
(counting their multiplicities) of a{—iz) on (0, 7V/2) is 
{~Y , so that in particular for odd £ the function a{~iz) 
must have at least one zero on the stated interval. Nu- 
merical work should urgently verify this last conclusion, 



as well as (|8.1| ), (8.2) in general, also in the context of the 
claim in [^Ipthat Rcif7' < N/2 for ^ = 2, but > N/2 for 
£ = 3. Further, it should be investigated whether WKB 
methods or those outlined at the end of Section VI can 



establish the sign of iq;(w) analytically also for ilu — + oo. 
In any case, in view of (5.3), in the latter limit those 
signs should be opposite for the RWE and ZE, unlike the 
situation for iuj I 0. 

The way in which throughout this article analytic 
continuations in r and in uj have turned out to be re- 
lated points to some interesting function theory in C^. 
While these pure-mathematical aspects have not been 
pursued, as their consequence and almost coincidentally 
we have arrived at a fairly complete theory of anoma- 
lous points and miracles in general, that is, independent 



of the specifics of the potentials (1.3), (1.4). Namely, in 
practice most potential tails will be either exponential. 



covered by the Frobenius method in Section II A, or al 



gebraic, covered by the two-step theorem of Section V A 



In particular, the latter, while comparatively intractable 
in other contexts because of their diverging Born series, 
have turned out not to lead to anomalous points. 

A similar degree of generality has not yet been achieved 
for describi ng m erging and canceling modes, as occurring 
in Section |VIl|, or in jl^ for the Poschl-Teller model. 



A model-independent study of merging modes is avail- 
able as the Jordan- block perturbation theory of , but 
there one starts by assuming the presence of a higher- 
order mode, which is subsequently split. In particular, it 
is not known whether the cancellation phenomenon is re- 
stricted to long-range potentials. The study of such issues 
will be facilitated by generalizing to long-range potentials 
the concepts of a QNM norm and more generally a gen- 
eralized inner product between two-component vectors 
(composed of both fields and their associated momenta) , 
which are known to be highly effective for finite-range 
potentials JlOf . 

The 'practical' significance of this work mainly lies in 
its results on the 'special' family of QNMs at small a, 
and in the guidance it gives to any subsequent numeri- 
cal investigations. Namely, there seems to be no reason 
why the rotation of an astrophysical black hole should 
ever vanish exactly. Furthermore, even for real-time nu- 
merical experiments one does not predict the 'special' ZE 
QNM to be independently observable: since the QNMs 



are not complete , the highly damped non-oscillating 
signal of the latter is in principle indistinguishable from 
the late-time tail caused by the branch cut in g [ p4| , 
consistent with the ZE's SUSY-equivalence to the RWE, 



which has been shown in Section VB to have no 'special' 
QNM. 

However, apart from the dependence of our small-a re- 
sults on the analysis at a = 0, there is no further jus- 
tification required to study the RWE and ZE closely, 
given that they have the same status in gravity as the 
Schrodinger and Dirac equations for the hydrogen atom 
have in quantum mechanics. Indeed, given its global sin- 
gularity structure, especially the RWE is the logical next 
step up in complexity from the hypergeometric equation 
encountered in the hydrogen problem |2^. In fact, a 
look at a figure of the Schwarzschild spectrum (e.g., [p^) 
shows that the special frequency is completely central, as 
the place from which the low- and high-damping branches 
somehow emerge. Hence, one has the hope that an im- 
proved understanding of O will lead to further insight 
into the spectrum in general, for instance by relating the 
two numerically observed branches to the two potential 
tails p8| , or more ambitiously, by providing accurate and 
general asymptotic formulae. 
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APPENDIX A: AN AFTERTHOUGHT ON a 

Our aquaintance with a has proceeded stepwise. De- 
fined originally as a set of proportionality constants in 
(^.2|) , immediately below at least locally on the NIA the 
function a{uj) was considered, focusing on the orders of 
its zeros. Subsequently, in Section VIII the function was 
studied globally on this axis, and simple ideas of real 
analysis were used to make a statement about the zeros 
of a. The logical elaboration is to continue a analyti- 
cally from the NIA, observing that (2.2) then yields the 
continuation of giiu) {gj-(uj)) to the fourth (third) quad- 
rant and beyond. One attractive property of a is that it 
is invariant under finitely supported perturbations of V. 
In fact this last result can be sharpened at once: a per- 
turbation ' 



cannot change a(w) for |Imaj| < A/2, 
and subsequently this invariance can be continued to the 
whole w-plane. 
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Let us presently derive a further important property: 
if the potential V{x) satisfies the decay property of Sec- 
tion V A with P — oo and if in addition it is single-valued 
-oo (so that V{x>xo) 



near x 



oo 

3=2' 



jX ^ ), then 



a{uj) is entire. Thus, for these potentials the definition 
( ^.2[ ) is even more fortunate than was apparent at first 
sight: from a two- variable multiple- valued function it cre- 
ates a one- variable non-branching function. 

For a proof, supplement g{x,uj) with a branch cut on 
the negative real a;-axis (in addition to the cut in uj on 
the NIA), set g{x,iv) = {g{x,uj),g{x,—uj))'^ , and define 
the monodromy map by 



g{-A-ie,uj) = M{uj)g{-A+ie,u!) 



(Al) 



(e,A > 0). Elementary properties are Aijk{i^) — 
M3-j,3-k{-io) U,k = 1,2) and M*iio) = M-\-lu*); 
this last relation involves a matrix inversion because for 
complex X one has g*{x,u!) = g{x*,—uj*), so that under 
complex conjugation the roles of —A ± ie are reversed. 



Furthermore, the two-step theorem of Section VA im- 
plies that A^22(w) = 1 in the fourth quadrant. To relate 
a to Al, observe that for a given ui, the two-step theorem 
yields the asymptotic form of g for a range of arg x obey- 
ing condition (C) of Section V A. Subsequently, g can be 
continued at fixed lu to the physical region x > by use 
of (Al) — effectively adding a third step to the two-step 



theorem — and comparison with (2.2) allows one to read 
off a. Thus, for instance for Rew < one finds 



(A2) 



making the relation between lu- (a) and x-continuations 
(Ai) maximally explicit [Q. Now let w = i\uj\, and 
consider the equality g{—A—ie,uj—j]) — A^ii(w— 77) x 
g{—A-\-ie,uj—r])-\-Mi2{tjJ—7])g{~A+ie, t]—uj) for 77^0. By 
the preceding one has A^ii(cj— 77) = 1, while in the limit 
g{—K—ie,uj) ~ g*{~K+ie,Lo) and g{—K-\-ie,rj—uj) 
g-[{—K-\-ie,—Lo). Since the latter function is ^ e^'^^ and 
thus real for ej,0, it f ollow s that Relim,,|o ■A4i2(t^— ??) = 
0. Combination with (A2) and the symmetry q;(— a;*) = 
—a*(Lo) then shows that lim^|o[Q^('^~?7) ~ Q^(^+'?)] — 0; 
so that a has no branch point at the origin. Since a 
Born calculation (cf. (B.l)) shows that a is also bounded 
near the origin, and since under the stated conditions on 
V the function cannot have singularities for finite oj (cf. 
Section VA), the derivation is complete. 

The above is a result of some power: while for algebraic 
potential tails the large- a; expansion of g furnished by the 
Born series diverges, one now sees that under the stated 
conditions on V the same series used as a small-o; expan- 
sion of a necessarily cotz verges. Also, in this situation 
one has explicit expressions for all branches of g; for in- 
stance, g[ijje~'^'^^) = [1 -I- a{uj)a{—uj)\g{ijj) -\- a{uj)g{—uj) 
if Rew > 0. Unfortunately, the only available exactly 
solvable example V{x>xo) = a{a'-{-l)x~'^ has a{uj) = 
2isin(7rtT), i.e., an w-independent constant which is not 



particularly instructive. Note that (pending a closer in- 
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